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Recent simulations and observations suggest that star clusters form via the as- 
sembling of smaller sub-clusters. Because of their short relaxation time, sub-clusters 
experience core collapse much earlier than virialized solo-clusters, which have simi- 
lar properties of the merger remnant of the assembling clusters. As a consequence it 
seems that the assembling clusters result in efficient multiple collisions of stars in the 
cluster core. We performed a series of A^-body simulations of ensemble and solitary 
clusters including stellar collisions and found that the efficiency of multiple collisions 
between stars are suppressed if sub-clusters assemble after they experience core col- 
lapse individually. In this case, sub-clusters form their own multiple collision stars 
which experienced a few colHsions, but they fail to collide with each other after their 
host sub-clusters assemble. The multiple collision stars scatter each other and escape, 
and furthermore the central density of the remnant clusters had already been depleted 
for the stars to experience more collisions. On the other hand, if sub-clusters assemble 
before they experience core collapse, the multiple collisions of stars proceed efficiently 
in the remnant cluster, and the collision products are more massive than virialized 
solo-clusters and comparable in mass to cold solo-clusters. 

Key words: galaxies: star clustersigeneral — methods: N-body simulations - galax- 
ies: star clusters: individual: R136 — open clusters and associations: individual: 
NGC3603, Westerlundl, Westerlund 2 



1 INTRODUCTION 

Young dense star clusters observed in the Milky Way 
and the Lar^e Magel l anic Cloud (L MC), e.g., R136 
Massev Hunt ei| 1 19981: ICrowther et aP 1 20H]| ) . NGC 3603 
Stolte et al I 200 61: 'Haravam a et al.l 20081) , "Westerlund 1 
■ Clark et al.l 200 5: Brandne r et al.l 20081: [Cennaro et al.l 
1201 ih and 2 ([Xscenso et al.1 l2007l : Rauw et al. 2007), are 
good samples for understanding the formation mechanism 
of dense star clusters. They are massive (- lO^Mo) and 
dense (> IO^Mqpc"^), and seem to be approaching (or 
might have experie nced) core collapse although they are 
young (< 4Myr) (Mack ev Gilmord ^3). For exam- 
ple, in R136 in the LMC, its high c ore density (> 5 x 
lO^Mopc-^) (|Mackev k Gilmordl 20031 ) and the existence of 
high-velocity stars (runaw ay stars ) escaping from the cluster 
(iBrandl et al . 2007; Eva ns et al.1 1 20 id : iBestenlehner et al.1 
I2OIII : Gvaramadze &: Gualandrid I2OIII) suggest that it ex- 
perienced core collapse (jFuiii Portegies ZwartI 1 20 111 ) . If 
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such a young massive cluster experiences core collapse, 
repeated collisions (so-called runaway collisions) of stars, 
and as a consequence the for mation of very massiv e stars 
(> lOOMcT)), are expected (jPortegies Zwart et "all Il999l : 
jPortegies Zwart fc^ McMillan I 2OO2I ). Such verv massive stars 
form through multiple stellar collisions could result in 
the formation of inte rmediate-mass black holes (IMBHs) 
(|Ebisuzaki et al.ll200ll ). 

The formation of IMBHs in dense star clusters via 
multi ple collisions has been stud i ed using A-body simula- 
tions (|Porteg . ies Zwart et al.ll 19991 . l2004l : ICiirkan et al]l2004l : 
iFreitag et al.l l20od ). and the results suggest that IMBHs 
with 10^ — 10^ M0 could be formed in such dense clusters. In- 
cluding stellar evolution, however, a high mass-loss rate due 
to the stellar win d of massive stars prevents the growth o f 
the massive stars (|Belkus et al.ll2007l : Iciebbeek et al.ll2009l ). 
A very high collision rate is required for such very massive 
stars to overcome the copious mas s- loss and never theless 
leads to the formation of an IMBH (jFuiii et al.ll2008l ). 

There are some mechanisms to enhance the growth rate 
of the very massive stars, but the most important factor is 
the moment of core collapse, tec- This short but high den- 
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sity phase is necessary for the cluster to become cohisionaUy 
dominated, which is critical for the collision rate of stars in 
the cluster. Earlier collapse times assist an efficient mass ac- 
cumulation because stars can start multiple collisions before 
the cluster starts to lose massive stars via stellar evolution. 

The core-collapse time is determined by the relax- 
ation time of the cluster if the cluster is initially virial- 
ized, and it is roughly 20% of the half-mass relaxation 
time, trh, with a Salpeter-type power- low mass functio n 
(jPortegies Zwart k McMillanI l2002l : iHeggie k HutI l2003l l. 
For most massive clusters, therefore, it is difficult to reach 
core collapse before the end of main-sequence lifetime of the 
most massive stars, which is ^ 3 Myr for > 40Mq. One way 
to achieve an early core collapse is kinematically cool ini- 
tial conditions. A sub-virial cluster ca n reach core collaps e 
faster than initially virialized clusters (lAllison et al.|[2009h , 
and that enables an efficient mass-growth via multiple col- 
lisions of stars. 

Another way to enhance the stellar-mass growth is by 
adopting mass segregated initial conditions. A mass func- 
tion causes massive stars to sink to the cluster center and as 
a result massive stars pile up in the cluster core. The mass- 
growth due to stellar collisions can be quite efficient, if mas- 
sive stars concentrate in the core. Initial mass segregation 
enhances the growth-rate o f the colliding stars (|Ardi et al.l 
I2OO8I : iGoswami et al.l I2OI2I I , and cool initial conditions re- 
sult in a high degre e of the mass segregation in a short time 
(| Allison et al.ll2009h . 

Another way for star clusters to reduce their 

core collapse time is by the as semblage of sub- 

clusters jAarse th Hillsl Il972l: [McM il lan etH] 
Moeckel Bonnelll l2009l: IYu et al.|[2Qnl: [Smith '^TaP 
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Fuiii et ahllioii ). The short relaxation time of sub-clusters 



compared to initially more massive clusters causes early 
mass segregation and core collapse. Since the memory of 
such earl y dynamical evo l ution is conserved in t he merger 
remnant (jMcMillan et al.1 l2007l : iFuiii et al.1 1201I hereafter 
Paper 1), the formation of star clusters by assembling them 
is also an effective way for efficient multiple collisions of stars 
in young star clusters. 

In paper 1, we found that 
of young dense star clusters 
clusters can successfully explain the mature character- 
istics of young massive star clusters such as R136 in 
30 Dor region. The age of R136 is only 2-3 Myr, 
but it shows dynamically mature characteristics, such 
as mass segregation, a high core densit y, and a wealth 
of high velocity es c aping stars (iMackev Gilmord l2003l; 
Brandl et aD l2007l: lEvans et al.1 boid : iBestenlehner et al.1 
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Gvaramadze Gualandris l201lh . However, the relax- 



ation time of R136 obtained from its current mass and radius 
is ^ 100 Myr (Mackey & Gilmore 2003), which is too long to 
have reached core collapse at its current age. In paper 1, we 
performed a series of A/'-body simulations of ensemble clus- 
ters and demonstrated that "ensemble" -cluster models can 
reproduce observations such as the core density, the fraction 
of high-velocity escapers, and the distribution of massive 
stars which experienced collisions, but "solo" -cluster models, 
which are initially spherical and virialized, fail to reproduce 
these observations. Furthermore, these characteristics of the 
ensemble models are also consistent with the characteristics 



of other massive young cluste rs like R136 in the L MC and 
NGC 3603 in the Milky Way (|Crowther et al.ll2010l ). 

If young dense clusters formed via assembling sub- 
clusters and have experienced core collapse, it is expected 
that repeating collisions can lead to the formation of 
very massive stars and possibly even IMBHs. In the ob- 
served young dense clusters, however, there is no evi- 
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In this paper, we perform a series of A/'-body simula- 
tions of solo and ensemble star clusters and demonstrate 
that the growth of very massive stars through multiple colli- 
sions is mediated by star cluster complexes. Our simulations 
show that the quick dynamical evolution of ensemble clusters 
does not always result in the formation of extremely- high- 
mass stars. When the assembling of clusters proceeds after 
each sub-cluster experiences core collapse ( "late- assembling" 
case), multiple-collision stars that form in each sub-cluster 
fail to coalesce to an extremely massive star, but leads to the 
formation of several very massive stars. Some of these very 
massive stars can escape from the cluster as high-velocity 
stars due to the three-body or binary-binary encounters. 
When the sub-clusters assemble before they experience core- 
collapse ( "early- assembling" case), the collision rate is en- 
hanced and the assembled cluster forms an extremely mas- 
sive star of - IOOOM0. 



2 METHOD AND INITIAL CONDITIONS 

We performed a series of A^-body simulations of solo clusters 
and ensemble clusters, that merge to a single cluster with a 
mass equal to solo clusters. For the ensem ble of sub-c lusters, 
we adopted two models. A: a King model (|Kindll966h with a 
dimensionless concentration parameter. Wo, of 2 and the to- 
tal mass Mci — 63OOM0, and B: a King model with Wo — 5 
and Mci = 2.5 x 1O^M0. The half-mass radii, rh, of these 
models are 0.092 and 0.22 pc, and the numbers of particles, 
N, are 2048 (2k) and 8192 (8k), respectively. The core den- 
sity is the same for both models (pc — 2 x 10^ Mq p c~^). We 
assumed a Salpeter initial mass function (IMF) (Salpeterj 
1955) between 1 and 100 Mq. We cah these models 2kw2 
and 8kw5. 

We distribute 4 or 8 of these sub-clusters in two dif- 
ferent initial configurations: spherical or filamentary. The 
former model stems from clumpy star formation in giant 
molecular clouds, and the latter is motivated by star for- 
mation in a filamentary gas distribution or shocked region 
of colliding gas in the spiral arms of a galactic disk. The 
clumpy star formation is init iated by observations of West- 



erlund 1 (lOennaro et al.ll201lh and R136 (ISabbi et al.ll2012h 
and simulations (jBonnell et al.ll201ll : ISaitoh et al.ll201lh . For 



the spherical models, we adopted 4 or 8 of models 2kw2 as 
sub-clusters, and distributed them randomly in a volume 
with a radius of rmax and with zero velocity. We varied rmax 
between 1 and 6 pc. For the filamentary models, we ini- 
tialized 8 individual 8kw5 model sub-clusters. We initialized 
these sub-clusters with two different initial mean separa- 
tions (models e8k8fl and e8k8f2), but with zero velocity. 
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Figure 1. Initial density profiles of single clusters 

The initial positions of the sub-clusters for these models are 
illustrated in Figure [21 All runs are summarized in table O 

For the solo models, we adopted two more initial con- 
ditions with Mci of 5.1 X IO^Mq and 2.0 x lO^M©. With 
the same mass function, these models have 16384 (16k) and 
65536 (64k) stars and are initialized using King models with 
Wo = 6 and 8, respectively. In order to obtain the same core 
density as that of sub-clusters, their half-mass radii are 0.32 
and 1.0 pc. We call these models as 16kw6 and 64kw8. In 
Table [T] we summarize the initial conditions, and we present 
their initial density profiles in Figure (3] We performed ad- 
ditional simulations of sub-virial (cold) initial conditions for 
16kw6, and an extra set of simulations in which we reduced 
the kinetic energy (velocity of each particle) to two-thirds 
and 10% of the virialized velocity. We call these models as 
sl6k-cool and sl6k-cold, respectively (see table [2]). 

The A/'-body simulations are performed using the sixth- 
order Hermite scheme with indiv idual timesteps with an ac - 
curacy parameter 77 =0.15-0.3 (" Nitadori Makinol l2008l V 
We adopted the accuracy parameter to balance speed and 
accuracy, and the energy error was < 0.1 % for all runs. Our 
code does not include special treatment for binaries, but the 
sixth-order Hermite scheme can handle hard binaries formed 
in our simulations (see section 2 in Paper 1). We took into 
account collisions of stars with a sticky-sphere approach and 
mass loss due to the stellar wind for stars with > lOOM© 
with a rate of 5.0 x lO-'^Moyr"^ (i Fuiii et al.|[2009l ). We ne- 
glected the mass-loss from stars with < lOOM© because it 
does not affect the results on the short timescale of our sim- 
ulations (< 5Myr). The stellar radii are taken from the zero - 
age main-sequence for solar metallicity (jHurlev et al.ll200Ql V 



3 SOLO-CLUSTER MODELS 

3.1 Virialized solo-cluster models 

We describe the results of the initially virialized solo-cluster 
models, which we will refer as the "standard" model. In 
Figure [3] we present the time evolution of the core density 
for models s2k, s8k, si 6k, and s64k. The core d e nsitie s are 
calculated using the method o f ICasertano HutI (|l985l ). We 
identify the moment when the cluster reaches the highest 
core density as the core-collapse time. The core-collapse time 




t (Myr) 

Figure 3. Time evolution of the core densities for solo clusters. 
The results are averaged in order to reduce the run-to-run varia- 
tions. 

measured from the simulations is tec = 0.29 ± 0.07, 0.71 ± 
0.11, 1.2 ±0.13, and 1.8 ±0.0 Myr for models s2k, s8k, sl6k, 
and s64k, respectively (see also Table [3}. The core-collapse 
time is cons istent with those obtained previous simulations 
(iGiirkan et al.ll2004r ). if we take into account the differences 
in t he mass range of the mass function. 

iGiirkan et al ] (|2004[ ) showed that the core-collapse tim e 
scales with the central relaxation time (|Heggie k Hutll2003l ): 

0.065cr^3D 

GHm)pAnA' 

Here C, (m), ctc, and, pc are the gravitational constant, the 
mean mass of stars, and the central velocity dispersion and 
density, respectively. Here In A is the Coulomb logarithm, 
and A is written as a function of the nu mber of particles 
asjjN and 7 :^ 0.1 for trh of star clusters (|Giersz Heggig 
I1994I I. We adopted the number of particles in the core as 
N. We find that the core-collapse time scales better using 
ou r definition than whe n we adopt 0.01 A", which is adopted 
bv IGiirkan et al ] (|2004l ). In our simulations, tcc/trc — 1 for 
models s8k, sl6k, and s64k, but tcc/trc — 0.5 for model 
s2k. For model s2k, however, trh is shorter than trc because 
the core radius exceeds the half-mass radius. If we adopt a 
shorter relaxation time, then tcc/trh ^ 1 for all the models. 

The core collapse of th e cluster initiates a cohision run- 
away in the cluster core (jPortegies Zwart et al.1 [19991 ). In 
Figure |4] we present the merger histories of the multiple- 
collision stars in the solo-cluster simulations s2k, s8k, si 6k, 
and s64k. In each model, one primary collision product 
(PCP) per cluster grows through repeated collisions of stars. 
In model s2k the mass-loss due to the stellar wind exceeds 
the mass-gain by the collisions, and therefore the PCP has 
lost all gained mass by the end of the simulation (5Myr). 
PCPs grow up to the maximum mass mmax ^ 4OOM0 via 
repeating collisions, but by the time it explodes the star is 
^ IOOM0. Here we define mmax as the maximum mass of a 
star reached during its lifetime as a result of collisions. 

PCPs are not the only stars that experienced collisions. 
In models s8k, si 6k, and s64k, we find secondary collision 
products (SCPs). In most cases SCPs experience only one 
collision (sometimes a few collisions) , but never grow as mas- 
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Table 1. Models of single clusters 



Model 


N 


Mel 


Wo 




Pc 


a 




trc 


Mcore/Mei 






(Mq) 




(pc) 


(M0PC-3) 


(km/s) 


(Myr) 


(Myr) 




2kw2 


2048 


6.3 X 103 


2 


0.097 


1.7 X 10^ 


11 


0.30 


0.58 


0.28 


8kw5 


8192 


2.5 X 10^ 


5 


0.22 


1.7 X 10^ 


15 


1.9 


0.92 


0.15 


16kw6 


16384 


5.1 X 10^ 


6 


0.32 


1.7 X 10^ 


17 


4.4 


1.1 


0.12 


64kw8 


65536 


2.0 X 10^ 


8 


1.0 


1.6 X 10^ 


19 


44 


1.8 


0.053 



a is the velocity dispersion. 



Table 2. Runs 



Model 




geometry 


(c^min) (pc) 


(sub-)cluster 




e2k4r3 


4 


spherical 


2.5 


2k2w 


3 


e2k4r6 


4 


spherical 


5.1 


2k2w 


1 


e2k8rl 


8 


spherical 


0.51 


2k2w 


2 


e2k8r3 


8 


spherical 


1.3 


2k2w 


1 


e2k8r5 


8 


spherical 


2.8 


2k2w 


2 


e2k8r6 


8 


spherical 


3.3 


2k2w 


2 


e8k8fl 


8 


filamentary 


2.8 


8kw5 


1 


e8k8f2 


8 


filamentary 


4.2 


8kw5 


1 


s2k 


1 






2kw3 


7 


s8k 


1 






8kw5 


6 


sl6k 


1 






16kw6 


6 


s64k 


1 






64kw8 


2 


sl6k-cool 


1 






16kw6 


2 


sl6k-cold 


1 






16kw6 


1 



The models are named according to the following rules; "e" and "s" indicate ensemble and solo models, respectively. For ensemble 
models, following numbers indicate the number of particles of sub-clusters and the number of sub-clusters. The last part indicates the 

initial configuration of sub-clusters; "r" and the following number mean spherical and the value of the maximum radius, rmax, "f" 
indicates filamentary initial configurations (see figure [2] for the initial positions of sub-clusters in these models). For solo models, the 
number indicates the number of particles. (cZmin) is the averaged distance to the nearest-neighbour sub-clusters, and A^run is the 
number of runs. sl6k-cool and sl6k-cold are the same model as sl6k, but the velocity of 67% and 10% of sl6k, respectively. 
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Figure 2. Initial position of ensemble models, e8k8fl (right) and e8k8f2 (left). We mimicked filamentary star forming regions. 
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Figure 5. Time evolution of the core density, pc, for models sl6k, 
sl6k-cool, and sl6k-cold. 



sive as PCPs, although SCPs sometimes exceed our adopted 
upper- limit to the IMF (lOOM©). The SCPs end up merging 
with PCPs (see bottom right panel in Figured]) or just lose 
their mass by stellar evolution (see top right panel in Figure 
This result agre es with previous numerical simulations 
(|Freitag et al.l l2006l ). 

We also find that the time when the PCPs reach their 
maximum mass TTlmax i ^max 5 is scaled by trc, and that 
tmax/trc = 2.3, 2.2, 2.2, and 2.6 for models s2k, s8k, sl6k, 
and s64k, respectively. 

3.2 Cold solo-cluster models 

Sub-virial (cold) initial conditions reach core collapse consid- 
erably earlier than virialized ones. Cold models have there- 
fore been suggested to explain the dynamic ally advanced 
appe arance of observed young star clusters ([Allison et al.l 
I2OO9I V In Figure [5] we present the core-density evolution of 
models sl6k-cool and sl6k-cold, which initially have 67% 
and 10% of the virialized temperature. These models reach 
core collapse much earlier than virialized models, and as a 
consequence multiple collisions start earlier and proceed at 
a higher collision rate. In figure [6] we present the mass evolu- 
tion of the PCPs for models sl6k-cold and sl6k-cool. Colder 
initial conditions result in a higher mmax of the PCPs. The 
high mmax is a result of the high collision rate, which is 
caused by the high density in the core (see Figure O . 

By the time the PCPs leave the main sequence (of 
^ 3Myr), their masses have been reduced considerably due 
to stellar mass-loss, which competes with the mass gain by 
collisions. In model sl6k-cold, the PCP grows quickly in the 
beginning of the simulation, but after 0.5 Myr the mass-loss 
rate due to the stellar wind becomes higher than the mass- 
growth by stellar collisions. In model sl6k-cool, the PCP 
stops growing at ^0.5 Myr because the mass- growth rate 
balances to the mass-loss rate, and then the PCP maintains 
its mass until the end of the simulation (3.5 Myr). Although 
the final masses of the PCPs are comparable in both sl6k- 
cool and sl6k-cold models, mmax of model sl6k-cold is twice 
as massive as that of model sl6k-cool. 

In Figure [71 we show mmax of the PCPs for all solo 
models. The maximum mass of the PCPs in models s8k. 
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Figure 6. Merger history of PCPs for models sl6k, sl6k-cool, 
and sl6k-cold. 



sl6k, and s64k is quite similar (^ 4OOM0) irrespective of 
Mci. One might expect that more massive clusters contain 
a larger number of massive stars and therefore a more mas- 
sive cluster can form a more massive PCP. In our simulation, 
however, the number of stars which merged into the PCPs 
and the mean mass of the merged stars are quite similar 
among these models (see Table [3|). By comparing models 
s8k, si 6k, and s64k, their (mcoi) and A^coi are quite similar 
even though their total cluster masses are different. If the 
collisions selectively occur among the most massive stars and 
the numbers of collisions are the same, larger clusters should 
have a larger mean collision mass (mcoi) because larger clus- 
ters contain more massive stars. However, the number of 
massive stars does not simply follow this relation. In figure 
[8] we plot the cumulative number distribution of massive 
stars with m > 5OM0 at the moment in which the mass of 
the PCP reaches mmax- The number of stars with > 5OM0 
within ^0.05 pc are similar (^ 20) among models s8k, sl6k, 
and s64k and slightly smaller for model s2k. In particular 
for the models sl6k and s64k, the distribution of massive 
stars preserves the initial distribution in the outer part of 
the cluster because the half-mass relaxation time exceeds 
tmax- The dynamical evolution in these models is driven on 
a timescale of trc, and they have similar core properties: trc, 
Mcore, and pc (see Table [1] and Figure [3]). 

In model s2k, on the other hand, trc trh and as a 
consequence the dynamical evolution proceeds throughout 
the entire cluster. Model s8k shows an evolution similar to 
that of model s2k: mmax/Md for model s8k is as high as 
that of model s2k. For these models trh 2 Myr, which is 
sufficiently short for massive stars in the outer part of the 
cluster to join the collisions in the core. Similar to the model 
s2k, models sl6k-cool and sl6k-cold can also gather massive 
stars from the entire cluster to the cluster center irrespective 
of their initial positions. In addition, these sub-virial models 
achieve very high density (see Figure [5]), which enhances 
the collision rate. The massive stars in model sl6k-cold are 
more concentrated towards the cluster center compared with 
model sl6k (see Figure [8]). 

Even though for model s2k the PCP can accumulate 
stars from the entire cluster population of massive stars. 
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Model s2k-6 Model s8k-1 




t (Myr) t (Myr) 

Figure 4. Mass evolution of PCPs (solid curve) and SCPs (dashed curve) for models s2k, s8k, sl6k, and 64k. Dotted line indicates the 
core-collapse time. Cross indicates the time when the star merged with more massive ones. 



their total number and mass still cannot compete with the 
population of massive stars in the more massive clusters. 
In these latter models, the maximum mass of the PCP is 
limited by the reservoir of massive stars, which manages to 
segregate to the core by the moment of the core collapse. A 
larger cluster mass therefore does not automatically lead to a 
massive PCP. As seen in Figures [4l and [6l the mass evolution 
of the PCPs in models s2k, s8k, and sl6k-cold show a clear 
peak in the middle of the simulation. In the later phase, 
when the collision rate decays, their mass-loss rate exceeds 
their mass-growth rate by stellar collisions. In models sl6k 
and s64k, on the other hand, they have not exhausted their 
reservoir of massive stars because their half-mass relaxation 
time is not shorter than the simulation time and therefore 
some of the massive stars still remain in the outer part of 
the clusters. 

We empirically obtained a relation that mmax = 
0.02Mci (dotted line in Figure [7|) for the low cluster-mass 
models (Md < 2 x 10^ M©) and the cold model. For mas- 
sive clusters, however, mmax is smaller than that accord- 
ing to this relation. For the most massive cluster (Md = 
2 X 10^ Mq), mmax is consistent with the result presented by 
IPortedes Zwart McMillanI (|2002l ). mmax = 0.002Mci. 
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Figure 7. The maximum mass of PCPs for solo models. Filled 
circles with error bars indicate models s2k, s8k, sl6k, and s64k 
from left to right. Cross and plus indicate sl6k-cool and sl6k-cold, 
respectively. Since the error bar for model sl6k-cool is smaller 
than the marker size, we do not plot the error bar. Dotted and 
dashed lines indicate ?Timax = 0.02A^d and ?77-max = 0.002Mci, 
respectively. 



4 ENSEMBLE-CLUSTER MODELS 

In section 3 we demonstrated that the results obtained from 
our solo-cluster models are consistent with previous numeri- 
cal studies. In this section we present the results of ensemble- 
cluster models, in which sub-clusters assemble to finally form 
one single cluster. In ensemble-cluster models, sub-clusters 



collapse on a timescale shorter than that for solo-clusters 
with the same total mass. Their further evolution is dom- 
inated by the dynamical evolution of the sub-clusters be- 
fore they merge. The conservation of the d ynamical states 
through the mergers (jVesDerini et al.ll2009l l drives the fur- 
ther evolution of the cluster merger products. As a result, 
ensemble clusters tend to experience core collapse consid- 
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Table 3. Summary of the results. 



Model 


mmax(M0) 


^max (Myr) 


Emerge (Myr) 


tec (Myr) 


^scPs(M0) 


(mcoiXMo) 


A^col 


e2k4r3-l 


287 


0.45 


0.2-0.87 


0.29 ±0.07 


375 


78.3 


8 


e2k4r3-2 


260 


0.78 


0.03-1.2 




454 


60.6 


15 


e2k4r3-3 


268 


0.61 


0.6-1.3 




139 


56.7 


12 


e2k4r6 


238 


0.57 


2.2-2.7 




743 


46.8 


16 


e2k8rl-l 


998 


0.80 


0.03-0.38 


0.29 ±0.07 





44.2 


45 


e2k8rl-2 


667 


1.35 


0.03-0.32 




160 


69.4 


22 


e2k8r3 


530 


0.86 


1.0 -0.75 




147 


80.2 


14 


e2k8r5-l 


334 


1.11 


0.47-1.9 




1192 


57.0 


25 


e2k8r5-2 


486 


0.78 


0.03-2.0 




651 


61.2 


19 


e2k8r6-l 


245 


0.59 


0.77-2.4 




1367 


51.0 


24 


e2k8r6-2 


274 


0.42 


0.03-> 3 




970 


45.5 


21 


e8k8fl 


1310 


1.40 


0.4-1.2 


0.71 ±0.11 


268 


73.0 


42 


e8k8f2 


659 


2.28 


0.8-2.0 




995 


88.4 


33 


s2k 


182 ±21 


1.3 ±0.6 




0.29 ±0.07 


16 ±40 


53.3 


4.6 


s8k 


399 ± 60 


2.2 ±0.2 




0.71 ±0.11 


149 ± 115 


63.2 


11.3 


sl6k 


431 ± 54 


2.6 ±0.9 




1.2 ±0.13 


54 ±77 


65.8 


13.2 


s64k 


488 ± 57 


4.4 ±0.2 




1.8 ±0.0 





66.0 


15.5 


sl6k-cold 


1064 


0.59 




< 0.02 





46.6 


40 


sl6k-cool 


707 ± 36 


2.55 ±0.75 




0.325 ±0.075 





51.1 


28.5 



T 




r (pc) 

Figure 8. Cumulative number distribution of stars with m > 
5OM0 at fmax for models s2k, s8k, sl6k, s64k, and sl6k-cold. 



erably faster than solo clusters which have initially similar 
properties to those of the merger remnant of ensemble clus- 
ters. In paper 1 we already showed that the quicker dynam- 
ical evolution of ensemble clusters can explain the mature 
characteristics of young dense clusters such as R136 and 
NGC 3603. Here we use that enhanced dynamical evolution 
to study the PCPs. The early dynamical evolution of ensem- 
ble clusters is similar to that of cold solo-clusters. One might 
expect that ensemble clusters also result in the formation of 
massive PCPs, but we will show that the early evolution of 
ensemble clusters is somehow more complicated. 

In Figure [9] we illustrate the schematic evolution of two 
typical evolutionary paths of ensemble clusters. We find that 



the most important parameter for the evolution of ensemble 
clusters is the moment of assembling, tens, compared to tec 
of sub-clusters. If tec > tens ("early assembling"), the PCPs 
in the remnant cluster grow efficiently by stellar collisions 
because the short relaxation time of the sub-clusters drives 
mass-segregation and core collapse faster than solo clusters. 
This evolution is similar to that of cold solo- clusters. 

If tec < tens ("late assembling"), each sub-cluster expe- 
riences core collapse before they assemble and form a PCP 
per individual sub-cluster. The mass of each PCP is lim- 
ited by the sub-cluster mass as we described in section 3. 
After the assembling of two or more sub-clusters, the PCPs 
formed in the sub-clusters sink to the center of the rem- 
nant cluster and interact each other. Most of them, how- 
ever, are scattered and ejected from the cluster because they 
tend to reside in hard binaries with a massive companion. 
The PCPs tend to be in the hardest binaries with the most 
massive stars when they formed in the sub-clusters. In each 
binary-binary encounter following a sub-cluster merger, two 
PCPs may collide although they are also ejected without ex- 
periencing a collision. Therefore, the majority of the PCP- 
binaries are scattered or ionized, and only one PCP-binary 
survives in the remnant cluster by the time the assembly is 
completed. The surviving PCP cannot continue to grow in 
mass because by that time the central density of the assem- 
bled clusters has been depleted due to the early dynamical 
evolution. 



4.1 Stellar collisions in ensemble clusters 

In Figures [10] and llll we present the mass evolution of 
PCPs and SCPs in ensemble clusters. The left and right pan- 
els show early and late assembling models, respectively. In 




Early assembling model (^cc > ^ens) 




Late assembling model (^cc < ^ens) 




Figure 9. Schematics picture of two typical assembling processes. Early assembling (tec > tens)' Sub-clusters assemble before they 
experience core collapse. The merger remnant is more mass-segregated than solo clusters which initially have similar properties to the 
merger remnant because sub-clusters have a shorter relaxation time than the solo cluster. After their assembling, the remnant cluster 
collapses and a massive PCP forms. Late assembling (tec < tens)- sub-clusters experience core collapse and form small PCPs before they 
assemble. After their assembling, however, the PCPs do not grow efficiently because most of them are scattered from the remnant cluster 
by binary-binary encounters. 



early assembling models, one massive PCP per remnant clus- 
ter grows after the assembling of sub-clusters. Even though 
some of the sub-clusters start forming PCPs before assem- 
bling, the PCPs merge after the host sub-clusters merged. 
In late assembling models, on the other hand, each sub- 
cluster grows its own PCP, but most of them do not collide 
with each other even after the assembling of their host sub- 
clusters. 

We find the reason for the difference between early and 
late assembling cases in the density evolution of these clus- 
ters. In figure [12] we show the time evolution of the maxi- 
mum number densities for ensemble and solo clusters. Here 
we plot the maximum value of the local density, which is 
calculated using six nearest neighbours. (Note that the max- 
imum local density does not trace the density of one individ- 
ual sub-cluster.) In early assembling cases, the density in- 
creases on the core-collapse timescale of the solo sub-cluster 
(model s2k) , but the maximum density is higher than that of 
model s2k and rather comparable to those of the cold models 
(models sl6k-cold, sl6k-cool). The evolution after the core 



collapse is similar to that of the cold models. The density 
gradually decreases and eventually becomes comparable to 
that of virialized solo-clusters (model sl6k). 

The density in late assembling cases also grows on the 
core-collapse timescale of the sub-clusters until a peak is 
reached at ^ 0.5 Myr. The density decreases as quickly as 
that of the solo sub-clusters (model s2k), which is different 
from early assembling cases. By the end of the simulations, 
the number density of the late assembling cases is an order 
of magnitude lower than in the early assembling cases. The 
relatively low density prevents the growth of PCPs in the 
late assembled clusters. The effect of the difference in the 
density can be seen in the number of stellar collisions, A^coi 
in TableO In early assembling models (e2k8rl-l and eSkSfl) 
and the cold solo model (sl6k-cold), A^coi = 42 di 2 and 
^max = 1100 =b 130, but in late assembling models (e2k8r5, 
e2k8r6, and e8k8f2) Ncoi = 24 ± 5 and mmax = 400 =b 150. 

In late assembling models (for e2k8r5 and e2k8r6), the 
maximum mass of the PCPs is 200-400 M©, but the mass 
of the PCP is similar to those of multiple SCPs, which 
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Figure 10. Top: Time evolution of the separation between sub-clusters projected onto a;-axis (full curves) and the collisions of PCPs 
(black dots) for models e2k8rl (left) and e2k8r6-l (right). The positions of the dots show the collision time and the the sub-cluster 
to which the star initially belongs. Bottom: Mass evolution of PCPs and SCPs for models e2k8rl (left) and e2k8r6-l (right). Crosses 
indicate the time when the SCPs merged to PCPs. Arrows indicate the time when sub-clusters merged. In all panels, the shaded region 
indicates the core-collapse time with error obtained from the simulations of isolated sub-clusters. 




Figure 11. Same as Figure [TO] but for models e8k8fl (left) and e8k8f2 (right). 
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Figure 12. Time evolution of maximum local number density 
(local densities of six nearest neighbours) for models sl6k, s2k, 
sl6k-cool, sl6k-cold, e2k8rl (early assembling), and e2k8r6 (late 
assembling). 



were PCPs in the sub-clusters. This feature is consistent 
with young dense clusters such as R1 36 in the LMC, which 
contains five > IQQMc:) mass stars (|Crowther et al.l I2OIOI : 
iBestenlehner et al]l201lh . although there is no evidence of 
any extremely massive stars with ~ lOOOM© . 



4.2 Maximum mass of PCPs in ensemble clusters 

As we show in section HTTl early assembling of sub-clusters 
results in the formation of a PCP, while late assembling 
forms a less massive PCP and multiple SCPs as massive 
as the PCP. In figure [13] we present the relation between 
^max/Mci and tenc/tcc of euscmblc models, where tec is the 
core-collapse time of the sub-clusters. Irrespective of the 
number of sub-clusters, the maximum mass of the PCPs 
decreases as the assembling time is delayed. 

In the left panel of Figure [Ml we show the relation be- 
tween TTimax of the PCPs and Md for both solo and ensemble 
clusters. (Note that for the solo clusters, the data is the same 
as that shown in Figure [7|). The PCP mass of early assem- 
bling models is higher than that of solo clusters with the 
same cluster mass and as massive as that of the cold model. 
In late assembling models, the PCPs is almost as massive as 
those of the solo clusters with the same mass. 

The difference in the maximum mass of PCPs is un- 
derstood if we take into account all the PCPs and SCPs 
in the cluster. In the right panel of Figure [M] we present 
the total mass of all the PCPs and SCPs in the cluster. 
The total masses are roughly located on the relation that 
mmax = 0.02Mci. This result suggests that the potential 
maximum mass of the PCPs is 2% of the cluster mass, al- 
though the value depends on the initial mass function and 
the mass- loss rate due to the stellar wind. The total mass of 
the SCPs is summarized in Table [3] as mscPs- These PCPs 
fail to merge with the most massive PCP and their mass 
will be lost from the cluster by escape or stellar evolution. 

In Figure 1151 we plot the radial distribution of PCPs 
and SCPs, which grows to > IOOMq. We combine the re- 
sults from several runs, separating them in the early and 
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0.00: 
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Figure 13. The maximum mass of the PCPs scaled by the total 
mass of the ensemble clusters as a function of the assembling time 
scaled by the core-collapse time of the sub-clusters for all e2k8 
and e2k4 models. Each horizontal line corresponds to one model. 
Vertical lines and crosses indicate the individual merger time of 
sub-clusters for four and eight sub-cluster models, respectively. 
The dotted line indicates tens /tec = 1- 



late assembling cases. While all the PCPs are located in the 
cluster core in the early assembling case, ^ 40% of the PCPs 
are ejected from the clusters or located in the outskirts of 
the cluster (> 10 pc) in the late assembling case. The num- 
bers of PCPs per cluster are on average 1.75 and 5.8 for the 
early and late assembling cases, respectively. In Figure [15] 
we also present the cumulative numb er distribution of stars 
with > lOOMcT) in the R 136 region (|Crowther et al.1 l20ld : 
IBestenlehner et al]l201lh . The number of such massive stars 
and their distribution imply that R136 experienced some 
late assembling, and observationally a sub-cluster has been 
found around R136 (|Sabbi et al.ll2012l ). 



5 SUMMARY AND DISCUSSION 

We performed a series of A^-body simulations of solo and en- 
semble star clusters and found that ensemble clusters evolve 
through typically two paths depending on their assembling 
time compared to the core-collapse time of the sub-clusters. 

In the early assembling case (tec > tens), the remnant 
clusters have dynamically mature characteristics (mass seg- 
regation and core collapse) compared to solo- clusters. The 
evolution of early assembling clusters is similar to that of 
sub-virial solo-clusters. The early assembling clusters expe- 
rience mass segregation and core collapse on the time scale 
of the sub-clusters, which is shorter than that of initially 
large solo clusters, and the short relaxation time of sub- 
clusters is conserved in the remnant clusters. This dynam- 
ically early evolution results in efficient multiple collisions 
of stars and helps the formation of extremely massive PCPs 
with - IOOOM0. 

In the late assembling case (tec < tens), the dynamically 
mature characteristics suppress the growth of massive stars 
via stellar collisions. In this case, the sub-clusters experi- 
ence core collapse individually and form their own PCPs, 
but the maximum mass of the PCPs in the sub-clusters is 
limited by the total mass of the sub-clusters. Even after the 
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Figure 14. The maximum mass of the PCP, mmax, in the cluster (left) and the total mass of mmax and the sum of the maximum mass 
of the SCPs, mscPs (right). 
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Figure 15. Cumulative distribution of PCPs and SCPs as a func- 
tion of the distance from the cluster center. Dashed and filled 
curves indicate early and late assembling models, respectively. For 
the early assembling models, we combined the data from e2k8rl- 
1, e2k8rl-2, e2k8r2, and eSkSsl (4 runs), and the average number 
of PCPs per run is 1.75. For the late assembling models, we com- 
bined the data from e2k8r5-l, e2k8r5-2, e2k8r6-l, e2k8r6-2, and 
e8k8s2 (5 runs), and the average number of PCPs is 5.8. Squares 
indicate the dis tribution of massive ( > lOOMcr)) stars observe d 
in R136 region (ICrowther et al.ll2010l : iBestenlehner et al.llioTlh . 
Since the observation is projected distance, we multiplied them 
byv^- For both the simulations and observations, we treat stars 
within 0.1 pc as at 0.1 pc because the distance is affected by the 
definition of the cluster center. 



sub-clusters assemble, the PCPs stop growing because the 
central density of the remnant cluster is already depleted 
due to the quick dynamical evolution of the sub-clusters. 
Since the PCPs in sub-clusters form massive binaries, they 
interact with each other in the remnant clusters. Some of 
them (SCPs) collide, but the others are scattered from the 
cluster by three-body or binary-binary encounters. In our 
simulations, 40% of the SCPs are ejected from the cluster 
or scattered to the outskirts of the remnant clusters. The 
SCPs sometimes escape with a high velocity (> 30km/s) 
and reach ~ 100 pc from the cluster within their life time 
(~ 3 Myr). The observed massive high- velocity stars such 
as VFTS 682 might be formed in this way (see also Paper 

1)- 



We also investigated the maximum mass of the PCPs 
and found that in ensemble clusters, the maximum mass de- 
pends on the assembling time of sub-clusters. In the early 
assembling models, the maximum mass of the PCPs is com- 
parable to that of sub-virial solo-clusters. In the late as- 
sembling models, however, the maximum mass is similar to 
that of the solo sub-clusters. The difference between them 
is mainly caused by the number of collisions. In the late as- 
sembling models, a larger number of SCPs are ejected from 
the cluster and fail to merge to the PCP than in the early 
assembling case. 



When the collisions of stars proceed most successfully 
(in early assembling and cold solo models), we find that 
the maximum masses of the PCPs reach ^2% of the total 
mass of the clusters even if we take into account the high 
mass-loss rate due to the stellar wind. Assuming an R136- 
like cluster of ^ 5 x lO'^Mo, the expected maximum mass is 
^ IOOOM0 . Such an efficient mass growth might result in the 
formation of IMBHs. For lower metalicity, the massive stars 
are p redicted to collapse directly to IMBHs (|Heger et al.l 



In late assembling cases, however, multiple smaller 
PCPs (IOO-4OOM0) are expected to exist inside or around 
the remnant clusters. These stars are in the mass range 
of type Ib/c supernovae (SNe) assuming solar metallicity 
(|Heger et al.| [2003V In recent observations of dense molecu- 
lar clouds in the central molecular zone in the Galactic cen- 
ter, several expanding shells were found, and the estimated 
total kinetic energy o f them is ^ 10^^ erg. (jTanaka et al.l 
I2OO7I : IOkaen il l2Q12l ). Especially, three major shells have 
a kinetic energy of ~ 10^^ erg, which corresponds to a hy- 
pernova explosion. A young dense massive clusters which 
is similar to our late-merger models might be embedded 
in this dense molecular cloud. Furthermore, escaping PCPs 
will explode up to ^ 100 pc from the host cluster. Ac- 
tually type Ib/c SNe associate with star forming regions 
(And erson et al ]|2010l : iLeloudas et°aDl201ll : lAnderson et al.l 
■2012 : ICrowthei]l2012l ). and for example Typ e Ic SN 2007gr 
is loc ated at ~ 7 pc from a young cluster ([Crockett et al.l 
I2OO8I ). 
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